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Abstract. We model the stochastic attenuation by Hi absorbers in the intergalactic medium 
fS| (IGM), such as LyQ Forest clouds, and absorbers associated with galaxies, such as Lyman 
— » Limit systems (LLS) and Damped Lyman Alpha absorbers (DLAs), and compute an ensemble of 
4 • 10^ attenuated Spectral Energy Distributions (SEDs) in the Johnson system for the spectrum 
■ of a galaxy with a constant star formation rate (CSFR). Using these, we asses the impact of 
the stochastic attenuation on the estimates of photometric redshifts for this type of galaxy by 
comparison with model SEDs that include only a mean attenuation. 
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. !• Motivation 

2^ The stochastic distribution of the (high column density) Hi absorbers causes a sig- 

' nificant scatter in the broadband colors of galaxies, especially at redshifts close to the 

\^ \ drop-out redshift in the corresponding passband where the photoelectric absorption is 

■ more severe, and this is expected to drastically affect the estimates of photometric red- 

\ shifts. 
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O ^ 2. Approach 

I We model the attenuation due to Hi along a random line of sight (LOS) using differ- 
ed • ential distribution functions constrained from observations l|Kim et al. 1997|l in a Monte 
^ Carlo fashion ( |Bershady et al. 1999| ) as described in Tepper Garcia & Fritze-v.A. (in 
prep.). We generate an ensemble of 4 • 10'^ lines of sight out to a given redshift Zem, 
each of them containing a random absorber population. For each LOS we calculate an 
^ . absorption mask, i.e. we compute the photoelectric and Lyman-Series line absorption 
■ ■ ' (as yet just for the first five Lyman transitions) caused by each absorber for a flat in- 
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put spectrum, modeling line profiles as in Tepper Garcia (2006) We compute model 



galaxy spectra corresponding to a CSFR for redshifts in the range 0.0 < Zem < 4.5 using 
the Evolutionary Synthesis code GALEV l|Bicker et al. 200411 . For a given redshift, we 
multiply each of the 4 • 10^ masks with our input spectrum and thus obtain an ensem- 
ble of an equal number of attenuated spectra. For each of these we compute a SED in 
the Johnson system. Using AnalySED IjAnders et al. 2006|l and a set of template SEDs 
that include only the mean attenuation IjJVladau 199511 for every redshift, we determine 
to which extent the redshifts of our simulated spectral energy distributions are recovered. 
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Figure 1. Estimate of the photometric redshift Zphot for an Sd-type spectrum with stochastic 
attenuation, using as input the same spectrum with mean attenuation. The sohd line shows the 
median, and the shaded area gives the 3cr-range around the mean for the entire ensemble of 
4000 SEDs. 

3. Results & Conclusions 

As can be seen in figure 1, we find that, using the mean attenuation only, the median 
photometric redshifts are severe underestimated by typically l^Zphot = 0.3, especially in 
the range z > 3.2. In addition, we observe a substantial scatter of Azphot ~ ±0.5. 

We predict that any method for determining photometric redshifts that only takes into 
account the mean attenuation at high redshifts should show a bias such as the one seen 
in Figure 13 We hence emphasise the need for an accurate modeling of the attenuation 
including stochastic effects, in order to correctly interpret in combination with evolution- 
ary synthesis models such as GALEV, the observations of high-redshift galaxies in deep 
multi-band imaging surveys. 
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